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Abstract. The quantum resonances occurring with 5— kicked atoms when the kicking 
period is an integer multiple of the half- Talbot time are analyzed in detail. Exact 
results about the momentum distribution at exact resonance are established, both in 
the case of totally coherent dynamics and in the case when decoherence is induced by 
Spontaneous Emission. A description of the dynamics when the kicking period is close 
to, but not exactly at resonance, is derived by means of a quasi-classical approximation 
where the detuning from exact resonance plays the role of the Planck constant. In 
this way scaling laws describing the shape of the resonant peaks are obtained. Such 
analytical results are supported by extensive numerical simulations, and explain some 
recent surprising experimental observations. 

PACS numbers: 05.45.Mt,03.65Yz,72.15Rn,42.50Vk 
1. Introduction. 

The present work is devoted to a detailed theoretical analysis of some of the quantum 
resonances occurring in the 5-kicked rotor, motivated by laboratory results obtained 
with experimental realizations of that system. 

1.1. Background. 

The (5-kicked rotor is a paradigm model in quantum chaology, and its theoretical 
importance is connected with the long-time properties of its evolution. Among these 
are effects of purely quantum origin such as dynamical localization ||, ||, ^, 0, §] 
and the quantum resonances @,|10|- Despite compelling numerical evidence, dynamical 
localization could not be mathematically proven until very recently||. The origin and 

I A proof for the case of sufficiently small kicking strength has been announced by J. Bourgain and 
S. Jitomirskaya while the present paper was being prepared for publication. 
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the main dynamical and spectral features of the quantum resonances are instead well 
understood in general mathematical terms [y, 0,[1^,[T^. They are observed when the 
kicking frequency is commensurate to the natural frequencies of the rotor and typically 
result in a quadratic growth of the energy asymptotically in time. 
Although the 5-kicked rotor is an abstract theoretical model, it has found an 
experimental realization using tools of atom optics, by a technique pioneered by 



Raizen and coworkers In present-day experiments laser-cooled Cesium atoms are 
periodically driven with a standing electromagnetic wave. As the frequency of the wave 
is slightly detuned from an internal atomic transition, a net force on the centres of 
mass of the atoms arises, proportional to the square of the driving field |1l4|,|T5i. The 
standing wave is periodic in space, so the atoms are subject to a periodic potential. 
The latter is turned on and off periodically in time, resulting in a sequence of short 
pulses, or "kicks" . In such experiments, the kicked atoms behave as individual particles 
with negligible interaction due to collision etc. |T6|, |1^ , so a single-particle theory is 
applicable. The nonzero duration of the kicking pulses sets a bound on the momentum 
range wherein the ideal (5-kicked rotor model is applicable. At large momenta the driving 
becomes adiabatic, leading to trivial localization in momentum for the classical and the 
quantum version of the rotor |jT8|,0- Important properties of the (5-kicked rotor have 



been experimentally reproduced: e.g., the exponential localization in momentum [|T3|,IT9 



(away from resonances). In addition the influence of noise and decoherence p0|, pl| , |22| , p^ , 
and the effect of gravity in the case when the kicking direction is vertical have been 
examined in experiments (with possible applications in atomic interferometry p^ , p5| , p6 
and quantum measurement theory P^, P^, PD|) . 

While dynamical localization is a robust phenomenon, the quantum resonances are 
rather sensitive instead, and the correspondence between theory and experiments is less 
appealing in their case. The following items represent incomplete or missing matches 
between the 5-kicked rotor theory and experimental results [^ ,|32|,^: 

(i) quantum resonances occur for the (5-kicked rotor whenever the kicking period r (in 
appropriate units) is a rational multiple of An, which is the Talbot time in the language of 
diffraction optics |3^,|3^,^. In experiments, only at the values r = 2n, An, 6tc, a special 
behaviour has been observed, in the form of narrow peaks in the energy absorption 
vs. r. This means that higher order resonances have not been resolvable within the 
experimental bounds. 

(ii) no quadratic energy growth at resonance has been experimentally observed. 
"Ballistic peaks" were reported by Raizen et al. [31], due to a tiny fraction of atoms 
which indeed show a quadratically increasing energy. This is in contrast to the bulk of 
the atomic cloud, which is instead frozen in a rather narrow distribution in momentum. 
The latter is, however, not exponential, at variance with the distribution found at non- 
resonant values of r, i.e. in the localized regime. 

(iii) an enhancement of the resonant peaks has been experimentally observed with 
increasing degree of decoherence. This is especially surprising, because decoherence 
is expected to drive the quantum system towards a more classical behaviour, hence to 
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inhibit purely quantum effects, such as the quantum resonances of the 5-kicked rotor. 
As to dynamical localization, its destruction due to decoherence was experimentally 
observed in ||,^,|21[|. 



The listed findings do not match with the theory of the 5-kicked rotor, and various 
physical reasons have been identified. The simplest, yet deepest of these is the starting 
point of the present paper: the experiments do not provide a realization of the 5-kicked 
rotor but of the 5-kicked particle instead, the difference being that while rotors move 
in circles, particles move in (approximately) straight lines. This difference is irrelevant 
on the classical level, however it has nontrivial quantum implications, that lie within 
the scope of the Bloch theory. The quantum dynamical localization is not impaired 
(localization lengths and fluctuation properties are affected, though) |]10|- Not so in 
what concerns the quantum resonances, that make the subject of this work. 

1.2. Outline, and statement of results. 

In this paper we analyze in mathematical detail the dynamics of atoms at quantum 
resonance and in its vicinity for the ideal case of 5-kicks. Our basic tool is the Bloch 
theory for quantum particles in periodic potentials. The basic constructions of that 
theory as applied to (5— kicked particles are shortly reviewed in section |^. Our main 
results are: 

(i) In the absence of decoherence, we show that the mean energy increases linearly with 
time (section |3.2.3| ) yet the motion is dynamically localized, in the sense that the 



momentum distribution settles in the course of time to a relatively narrow steady 
state distribution (section |3.2.2 ). This seeming contradiction is resolved by the slow 



algebraic decay of the tails of the steady state distribution. We give an integral 
expression (eq. (iOD,(|ll),([l3)) of the latter distribut ion, whence we derive exact 
bounds (^) and large-momentum asymptotics (^2]) . 

(ii) We derive a scaling law (section ^ that describes the shape of the resonant peaks 
which are observed in the energy vs. kicking period curves. The scaling law (^) 
is then demonstrated by numerical simulations. It is derived by means of a quasi- 
classical approximation (subsection [4.1|) , where the role of the Planck constant 
is played by the detuning from exact resonance. This technique was introduced 
in [p5| , |26[| . It allows to describe the near-to- resonant quantum motion in terms of a 
Standard Map, which is different in parameter values from the one that is obtained 
in the usual classical limit of the (5-kicked rotor. The quantum resonance exactly 
corresponds to the main c/asszca/ resonance of this standard map. The stable island 
associated with the latter resonance accounts for the structure of the resonant peaks. 
This analysis closely links quantum resonances with classical nonlinear ones, thus 
closing a gap between different meanings of the term "resonance", which were 
deemed to be unrelated. 

§ Similarly noise due to Spontaneous Emission affects non-dispersive quantum wave pacfcets formed in 
liighly excited Rydberg atoms p7| , |38|] . 
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(iii) We devise a model of stochastic (5-kicked rotor dynamics, where the external 
noise is modelled after the Spontaneous Emission (SE) effects that are used 
to induce controlled decoherence in experiments (section |^). Our statistical 
assumptions (subsection |5.2| ) are a compromise between physical adherence and 
mathematical convenience, and allow to derive simple exact results, that expose 
the key mechanism whereby the quantum resonances are affected. We explicitely 
compute the growth of the mean energy in time (eq. (|6^) ). We prove that no 
steady state momentum distribution is any more attained, and that, in contrast, 
the momentum distribution evolves into a Gaussian distribution which spreads 
diffusively (subsection |5.2.4|) . By checking such results against realistic numerical 
simulations we demonstrate that their validity does not crucially depend on our 
technical assumptions (subsection |5.2.5| ). So they are relevant for a wide range of 
experimental situations. 

(iv) The stochastic 5-kicked rotor dynamics provides a remarkable example where the 
classical and quantum effects of the coupling to the environment can be separated. 
This is achieved by a stochastic gauge transformation (section 

(v) We provide a description of the nearly resonant dynamics in the presence 
of decoherence (section |5.3|) , by means of the small- detuning quasi-classical 
asymptotics (point (ii) above). In this way we derive a scaling law for the structure 
of the resonant peaks in the presence of decoherence (eq. (|79|) ), which is then 
demonstrated by numerical simulations. 

(vi) We explain the experimentally observed enhancement of the resonant peaks in the 
presence of SE (section ^). At a low degree of decoherence, the linear growth 
of energy at resonance (eq. (0)) is but slightly changed with respect to the 
decoherence-free case, so decoherence cannot by itself explain the much stronger 
enhancement of the resonant peaks that was experimentally observed. In order to 
understand the latter we analyse differences which experiments inevitably impose 
with respect to the ideal models. Experimental cutoffs in momentum exist, due on 
one hand to the finite duration of the pulses and on the other hand to the signal-to- 
noise ratio. Using the results of point (i) we argue, and numerically demonstrate, 
that the latter type of cutoff is responsible for the observed enhancement. The 
coherent energy growth at resonance is due only to the fastest atoms, hence it is 
stopped as soon as such atoms escape the finite bounds of experimental observation. 
This peculiarity is rapidly destroyed by decoherence, which calls for all the atoms 
to participate in the mean energy growth. 



2. Kicked atoms vs. Kicked Rotors. 



We consider a one- dimensional model for a kicked atom of mass M, subject to time- and 
space-periodic kicks, with periods T and vr/c^^, respectively, is the wave number of 
the kicking field whose maximum strength is given by k. We rescale momentum in units 
of 2hkL, position in units of (2/ci)~^, mass in units of M. Energy is then given in units 
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of ^^(2fci)^/M, time in units of M//i(2fci)^, and the reduced Planck's constant equals 
1. The dynamics of the kicked atoms is induced by the following Hamiltonian p6|,p5[1, 
which depends on the continuous time parameter t': 

H{t') = — + kcosiX) J2 5{t'-tT), (1) 

t=~oo 

where X,P are the position and the momentum operator respectively, k = k/H, and 
the kicking period r = hT{2kiY / M . The state evolution of the atom from one kick to 
immediately after the next kick is determined by the unitary Floquet operator: 

U = e-ifccos(X)g-irPV2_ (2) 

Iterated application of the one-cycle operator U yields the dynamics of the atom in the 
discrete time given by the kick counter, which we denote by the integer t. We further 
denote \ip) the state vector of the atom, and ipi^) = {^14'), i'ip) = {p\4') the wave 
functions in the position and in the momentum representation, respectively. 
The model (0) differs from the kicked-rotor model because the particle does not move on 
a circle but on a line instead. A link between the two models is generated by the spatial 
periodicity of the kicking potential: the Floquet operator (0) commutes with spatial 
translations by multiples of 27r. As is well-known from the Bloch theory, this enforces 
conservation of the Quasi- Momentum (QM). In our units, QM is given by the fractional 
part {p} of the momentum p and will be denoted (3, [0 < (3 < 1). For a sharply defined 
value of quasi-momentum, the wave function of the atom is a Bloch wave, of the form 
exp{iPx)ipi3{x), with '?/'/3(x) a 27r— periodic function. The general particle wave packet is 
obtained by superposing Bloch waves parametrized by the continuous variable P G [0, 1): 

ij{x)= f d(3e'P^M^). (3) 

Denoting 6 = x mod(27r), we in turn have 

M0) = ^J2'^in + f3)e'-'. (4) 
yZn — 

* n 

In the special case when the state of the particle is a plane wave with momentum po, 
then 



ijjf.ie) = {V27r)-'6{(3-(3o)exp{znoe) , (5) 

where /3o = {po} and uq = [po] are the fractional and integer part of po, respectively. 

For any given (3, ippiO) may be thought of as the wave function of a rotor with 
angular coordinate 9, henceforth called /3— rotor. We denote the corresponding state of 
the rotor by iV'/j). From (0) and (|), it follows that, while \ip) evolves into U lip), li'/s) 
evolves into Ujs {ipp), with 

Up = e-*fc-W e-*i(^+«', (6) 

where A/ is the angular momentum operator: in the 6'— representation. A/ = —id/dO 
with periodic boundary conditions. The Floquet operator differs from the Floquet 
operator of the standard 5-kicked rotor by the phase 13. In previous studies of this rotor 
variant, (3 was typically regarded as an external Aharonov-Bohm flux [09|] . 
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3. Quantum resonances. 

The time evolution determined by the Floquet operator (^) preserves many of the 
dynamical properties of the standard (5-kicked rotor (corresponding to /5 = 0), but 
not necessarily the quantum resonances. These only occur for special values of quasi- 
momentum |T^. Whenever r = Anp/q {p,q mutually prime integers), and /3 = m/2p 



with m an integer such that < m < 2p, the operator (^) commutes with translations in 
momentum space by multiples of q, leading to bands in its quasi-energy spectrum. For 
special values of q, the bands may be flat, i.e., of zero width; such is the case, e.g., for 
g = 2, m an even integer. In typical cases, however, the bands are not flat and result in 
absolutely continuous quasi-energy spectra, that enforce ballistic growth of momentum. 
The width of the bands rapidly decreases as the order q of the resonance increases 0,[1O 



so ballistic motion is observable only after quite long times. This makes higher-order 
resonances experimentally hardly detectable, and provides an already well-known answer 
to the first problem (i) stated in the Introduction (subsection |1 . 1|) . Further reasons 
hindering experimental observation of higher resonances are pointed out in section ^. 

3.1. (3— rotor dynamics. 

We focus on the main resonances g = 1,2 in the following, i.e., we set r = 27r£, with 
I a positive integer. Inserting this value for the kicking period into (|^), and using the 
identity exp(— ivm^^) = exp(— ivmf), we obtain (apart from an n- independent phase 
factor) 

where ^ = ttI{2(3 ± l)mod(27r) will be taken in [— tt, tt). The 2nd operator on the rhs of 
(0) will be denoted 7^(0- the 6'— representation it acts according to: 

{nmpm=ijp{e-i). (8) 

The state after the t— th kick is then given by 

(W*^^)(^) = e-'=^(^'«'%(^^-tO, (9) 

with 

i-l 

F{e, ^,t) = Y, cos{e -sO = \Wt\ cosie + ^igiWt)) , (lo) 

s=0 

where Wt = Wt{^) := X]l=o^ denote by n the eigenvalues of the angular 

momentum Af. Then, in the J\f representation, the state (|^) reads (after changing 
variable from 6 to 6+aig{Wt)): 

(n |W* = e*"^^g(^'^ r"^e-^"l^*l^°^(^V/3(^-te-arg(iyi)) .(11) 



^27r 

If the initial state of the particle is a plane wave (|^) of momentum po = uq + i3o, then 
^ takes the constant value = 7r^(2/3o — 1)- Substituting in (|TT]), and computing the 
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integral by means of formula ( [li).l| ), the momentum distribution for the /^q— rotor at time 
t is: 



where J„ 
identity 
by 2): 



PMno,Po) = Jlnoi^m), (12) 

is the Bessel function of first kind and order n. Using the Bessel function 
, one computes the expectation value of p"^ (or the energy by dividing 



(no + /3o)' + ^fcWt 



Explicit computation of (|TOD yields: 

sin(teo/2 



if eo 7^ . 



(13) 



(14) 



sin(eo/2) 

If ^0 = then \Wt\ = t. Then the distribution (|12|) spreads linearly in time, and the 
average kinetic energy increases like k'^t^/A. The corresponding spread over the discrete 
momentum ladder is the same as for a free particle under the time evolution generated 
by the discrete Laplacian operator, = corresponds to Pq = 1/2 + n/i mod(l), 
n = 0,1, ..,£ — 1, the resonant values of /3 which were mentioned in the beginning of 
the present Section. For any other value of (3, the distribution changes in time in a 
quasi-periodic manner. It oscillates in time with the approximate period tc^q^, inverse 
to the detuning of (3q from the nearest resonant value. At any time t, it is negligible at 
\n — no\ > k\ csc(^o/2)|. 



3.2. Incoherent ensemble of atoms. 

If the initial state of the particle is a wave packet, then it is a coherent superposition of 
continuously many plane waves with different quasi-momenta, which are non-resonant 
except for a finite set of values. It can then be proven that the asymptotic growth of 
energy in time is proportional to k'^t/A ||2^. Here we consider the case when the initial 
atomic ensemble is an incoherent mixture of plane waves. Numerical simulations based 
on such choices of an initial state have shown satisfactory agreement with experimental 



data ||0,^. The initial momentum distribution shall be described by a density f{p)- 
We can equivalently consider an ensemble of /9— rotors with /5 distributed in [0, 1] with 
the density: 

+ 00 

/o(/3)= + (15) 

n=~oo 

In the case when f{p) is Gaussian with standard deviation cr (a reasonable assumption 



for the Oxford experiments according to |33|), /o(/5) is a Theta-function, and 
Poisson's summation formula yields: 



/o(/3) 



1 + 2e 



cos(27r/5) + 0{e- 



(16) 



For a > 1, that is relevant for the present day experiments, it is practically 
indistinguishable from the uniform distribution fo{f5) = 1. Each /5— rotor is described 
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by a statistical state, which attaches the probabihty fo{l3)~^f{n + j3) to the momentum 
eigenstate |n). The momentum distribution P{p,t) of the particle at time t is obtained 
as follows. For any given G [0, 1), averaging (p!2[) over the different no of the initial 
distribution yields the momentum distribution P{n,t\j3o) for the /3o— rotor. Weighted 
by /o(/9o); this is the same as the momentum distribution P{p, t) of the particle over the 
ladder p = n + Pq (/3o fixed, n variable). The on-ladder distributions corresponding to 
different /Sq combine like a jig-saw puzzle in building the global momentum distribution 
for the particle. The result is a complicated function of p that oscillates on the scale 1/t 
(see the dependence of (p!^),(|l4|)). Nevertheless, in the average, this distribution 
evolves into a steady state distribution. This may be shown either by time-averaging, 
or by coarse-graining, as done in the following subsections. 



3.2.1. Time- averaged distribution. As t ^ oo, the time average of the distribution 
(|l^) tends to whenever f3o has a resonant value, because the ballistic flight of atoms 
with a resonant quasi-momentum results in vanishing probability of being found in any 
finite momentum interval. For any non-resonant value of jSo, the distribution (^) tends 
in time-average to a distribution P*{n\no, Po). In fact, for non-resonant /3o the time 
average of (|12D is just the average of J^_^^{k sm{a) csc(^/2)) along the trajectory of the 
shift a — s> a + ^o/2(mod27r), starting from a = 0. If (3o is a non-resonant rational, then 
the trajectory is periodic and the time average trivially converges. If /3o is irrational, the 
shift is ergodic, so the time average converges to the phase average. The time-averaged 
momentum distribution for the /3o-rotor is thus given, for irrational (3q, by the following 
integral expression: 

1 f'^'^ 1 
P*{n\no,Po) = — J da J^_^^{ksm{a)csc{7ii{po- -)))■ (17) 

It follows that the momentum distribution P{p,t) for the atom converges, in time av- 
erage, to a steady-state distribution P*{p). This is obtained by averaging over no the 
time-averaged distribution P*(n|no,/?o) discussed above for the /5— rotors, and by com- 
bining the on-ladder distributions P*(n|/5o) thus obtained. Note that the rhs of ([T7|) , 
albeit a continuous function of (3q, does not necessarily coincide with the Ihs when [3q 
is rational. Hence P*{p) has, strictly speaking, a dense set of discontinuities formed by 
values of p with a rational non-resonant fractional part Nevertheless (|l^ is adequate 
for computing averages of smooth functions of p. A time-averaged distribution (p!7| ) is 
shown in Fig. |l| for ^ = 1 and for the case of an initially fiat momentum distribution 
in [0, 1). In this case uq = 0, and the integral ( [T7| ) was numerically computed for each 
n, Pq. The distribution is not seen to vanish at n + 1/2 because the computational grid 
in p was chosen at random, thus avoiding simple rational numbers. Discontinuities at 
integer p are clearly visible. They are due to the discontinuity of the initial momentum 
(box) distribution at p = 0, 1, and are not related to the above discussed ones. In the 
momentum intervals n < p < n + 1 with n ^ k the distribution P*{p) is negligibly 
small except in a neighbourhood of width 2/n around n + 1/2. This is due to the 
faster-than-exponential decay of the Bessel functions when their order increases beyond 
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the argument. At p = n + 1/2 the distribution is exactly zero. Thus in all intervals 
n < p < n + 1, P*{p) has the form of two narrow peaks, situated symmetrically wrt 
n+ 1/2. As shown in the next subsection, the total probability in such intervals decays 
like 1/n^, so the height of the twin peaks decays like 1/n, as is indeed clear from Fig. 
Cases with i > 1 may be discussed along similar though not identical lines. 



3.2.2. Coarse-grained distribution. The second way of removing the fast oscillations is 
replacing P{p,t) in each interval n < p < n + 1 by its integral Pn{t) over that interval. 
This corresponds to using a bin size 2hkL for the observed distributions. Assuming f{p) 
to be coarse-grained itself, the new distribution is approximately computed in the form: 

P„(t) = ^M„_„(t)/(m), (18) 

m 

where 

r dr 

M,,{t)= j dpjl{k\Wt\) = J —J^{ksm{tx)csc{x)) 

27r 

— ^2_,— (A;sin(x)csc(xt-^ + 27rrt-^)) . (19) 

-■^ r=0 ^ 

In the limit when t ^ oo and 27rr/t a, the sum over r approximates the integral over 
a, and (|19D converges to 

dx I da Jl{ksm{x)csc{a)) . (20) 

-TT Jo 



The steady-state coarse-grained distribution P* is then obtained by replacing (|20D in 
(p!8|). Note that it is the coarse-grained distribution itself, and not just its time-average, 
that converges to the steady-state distribution. We do not know if the double integral 
may be computed in closed form. In [Appendix A| we prove the following (non-optimal) 
estimate, valid for N > k: 



EK.2g ""A^i^(24 . (21) 

\n\>N ^ ^ ^ ^ 

Using this estimate it is easy to compute that for k > 1 the total probability carried 
by states \n\ > Ak is not larger than 0.31 . . .; so, at large k, the distribution is rather 
narrow as compared to the exponentially localized distribution which is observed far 
from resonance, because the width of the latter scales like fc^. In [Appendix A| we further 



prove that the distribution ([20|) has the following large- |n| asymptotics: 



Ak 

M* — as |n| ^ oo . (22) 

Such a decay carries over to the coarse-grained distribution P* whenever f{p) is fast 
decaying (e.g., like a Gaussian). 

The convergence of the coarse-grained distributions to the steady-state distribution 
is illustrated in Fig. |], that were produced by numerically simulating the evolution 



QR 



10 



of statistical ensembles of particles, with an initial Gaussian momentum distribution. 
The central part of the distribution quite early stabilizes in the final form of a narrow 
peak of width ~ k. Away from this peak, the algebraic tail oc develops over 
larger and larger momentum ranges as time increases in the wake of two symmetric, 
tiny "ballistic peaks"!], that move away linearly in time. The fall of the distribution 
is quite steep past such peaks. This is easily understood from the first equation in 
(|19|): since | sin(ta;) csc(a;)| < 7rt/2, at > 7ckt/2 the integrand decays faster than 
exponentially. The distribution in Fig. ^ (b) has stabilized to the limit distribution over 
a broad momentum range. Apart from the far tail, where the moving peak structure is 
still apparent, the distribution follows the asymptotic decay (|2^) already for \n\ > 15. 
Hence, using (^), the total probability on states \n\ > 40 is ~ 8 x 10~'^. 
From the above analysis, we obtain that all moments oip of order > 1 diverge ast ^ oo, 
in spite of the onset of the stationary distribution, due to the slow algebraic decay of 
the latter. For the case of the 2nd moment, which is (apart from a factor 2) just the 
mean energy of the ensemble, the growth is actually linear in time, as we shall presently 
show. 



3.2.3. Average kinetic energy. The mean kinetic energy of an ensemble of rotors at 
time t is obtained by averaging (|13|) over the initial momentum distribution: 

P /•^,,,,,, sin^(t7r£(/?-l/2)) 
4 Jo sm (7r£(/3 - 1/2)) 

As t — s> oo, the fraction in the integrand, multiplied by C./t, tends to a periodic 6 function 
of {(3 — 1/2) with period l/C. Thus (^) has the following t ^ oo asymptotics: 



S{E{t)] ~ ^{^(0)} + y E /o(/?^^'^ (24) 

where (3^^^ = 1/2 + j/i mod(l) are the resonant quasi-momenta. In the case when 
/o(/?) = 1, i.e., it is uniform in [0, 1), this formula is exact for all times t (by ( |E.7| )). It is 
practically exact at all times for an initial Gaussian distribution with width a > 1 around 
the origin; for, indeed, fo{(3) is nearly uniform in that noted in section 13.21. 

Higher order energy moments may be likewise computed. For large, finite t the 
decay of the distribution is truncated ai n kt and is thereafter replaced by a much 
faster decay. Consequently, the variance of energy increases like t^. The increase of 
other moments may also be estimated in this way. 



4. Dynamics near to Resonance. 

The experimental observation of quantum resonances is mainly based on measuring 
the average energy of the kicked atoms after a fixed time tabs- For all sufficiently 

II This denotation is borrowed from |^l| , where experimental observation of such structures was 
reported. 
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irrational values of the kicking period r, the theory of the ideal kicked rotor model 
predicts localization, that is, on increasing tobs beyond a break-time t*, the observed 
energy values should not increase any more. In contrast, resonant values of r lead 
to unbounded growth of energy^f] with tobs, as analyzed in previous sections for the 
particular cases r = 27r£. Hence, if tobs is significantly larger than t*, then a scan of 
the measured energy versus the kicking period r yields plots similar to the (numerically 
obtained) one shown in Fig. ^(a), where peaks are clearly observed at the resonant values 
r = 2TTi. For obvious continuity reasons, such peaks have a width, determined by the 
finite value of tobs- In the ideal case, they would shrink on increasing tobs, and further, 
narrower peaks associated with higher-order resonances would appear. In this section 
we derive a description of the structure of the peaks around r = 27r£, based on a finite 



time, small-e asymptotics, where r = 27ri + e. This technique was introduced in ||25|,[26 
We in particular find that the width of the resonant peak scales like (ktlf^^)'^, so that 
at large tobs the peak is much narrower than the naive expectation oc I /tobs- 

4.1. e— quasi classical asymptotics near resonance. 
We rescale k = k/\e\, and define 

i = \e\Af = -i\e\^ , 7i^(/,t) = ^sign(e)/2 + /(7r£ + r/5). (25) 
Then the Floquet operator for the /5— rotor may be rewritten in the form: 

W^(t) = e-R^^°^('~) e"^^^ (26) 

If |e| is regarded as the Planck's constant, then (p5D,(p6D is the formal quantization of 
either of the following classical maps: 

It+i = It + ~k sm{et+i) , 9t+i =9t±It + ni + Tp mod(27r) (27) 

where ± has to be chosen according to the sign of e. We stress that "classical" here is 
not related to the h ^ limit but to the limit e ^ instead. The small— |e| asymptotics 
of the quantum /3— rotor is thus equivalent to a quasi-classical approximation based on 
the "classical" dynamics (0), that will be termed e-classical in the following. Changing 
variables to J = ±7 + 7i£ + r/5, = 9 + n{l — sgn(e))/2 turns the maps ( pT] ) into a single 
Standard Map, independent of the value of /3: 

Jt+i = Jt + ~ksmi^t+i) , A+i='&t + Jt- (28) 

This will be called the e-classical Standard Map (eSM) in what follows. In Fig. § 
quantum energy curves vs. r in a neighbourhood of r = 27r are compared with energy 
curves computed using the e— classical map (|27|) . For any given particle in the initial 
ensemble, the map (pTj) with j3 equal to the quasi-momentum of the particle was used 
to compute a set of trajectories started at I = no\e\ with homogeneously distributed 
9 G [0,27r). The final energies e~'^If/2 at t = tobs of the individual trajectories were 

^ Unbounded growth was also proven for a dense set of close-to-commensurate values of r 
Extremely long times are however required to resolve such arithmetic subtleties. 
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averaged over 6, f3, tiq with the appropriate weights. This is equivalent to using the 
eSM in all cases, with different initial ensembles Jq = const = ±r2o|e| + Tci + rfS^. As 
(5q is varied, such ensembles sweep the full unit cell of the eSM, so sampling different 
(3qS amounts to probing different regions of the e— classical phase space. The average 
energy {Et) = e~^{If)/2 is plotted vs. r = 27r + e in Fig. ||, along with results of the 
corresponding quantal computations. The main qualitative features emerging of Fig. ^ 
are: (i) On a gross scale the curves are shaped in the form of a basin with a high, narrow 
spike in the centre, closely flanked by a much smaller peak on either side, (ii) Quantum 
and e— classical curves nicely agree at small |e|, in particular the structure of the spike 
is the same. Their behaviour at large |e| is qualitatively similar but quantitatively 
different. 

This overall qualitative behaviour may be explained in e-classical terms, and an 
approximate scaling law for the t, k, e dependence of the mean energy close to resonance 
may be obtained, as shown in the next section. 



4- 1.1. e- quasi- classical analysis of the resonant peaks. The e— classical standard map 
is different from the map obtained in the classical limit proper /i — of the kicked 
rotor. In particular, if 27rik > 1, then the classical and the e— classical dynamics are 
at sharp variance whenever k < 1. In the former unbounded diffusion occurs; in the 
latter the dynamics is quasi-integrable instead, and the e— classical trajectories remain 
trapped forever in-between KAM curves. It is exactly the deep changes which occur 
in the e— classical phase space as r is varied at constant k that account for the energy 
vs. r dependence at fixed time t = tabs- In the following discussion we assume for 
simplicity an initially flat distribution of p G [0, 1]; then Jq = 0, and Jq = ni + t/Sq with 
Po uniformly distributed in [0,1). Without loss of generality we also consider i = 1. 
Hence if |e| <^ 1 then Jq is uniformly distributed over one period (in action) (tt, Stt) of 
the eSM. 

Since Jj = ±/t + ir + r/3, and Jq = 0, the mean energy of the rotor at time t is: 

(Et,,) = e-^I^)/2 = e-'{{6JtY)/2 , 6Jt = Jt - Jo ■ 

The exact quantum resonance e = corresponds to the integrable limit of the eSM, 
where 6Jt = 0. However, (Ef^^) is scaled by e~^, so in order to compute it at e = one 
has to compute 6Jt at first order in e. This is done by substituting the 0-th e-order of 

the 2nd equation into the 1st equation of (128|) . This leads to: 

t-i 

6Jt = |e|A;^sin(0o + ^os)+r(e,t) (29) 

s=0 

where r(e, t) = o(e) as e ^ at any fixed t. The energy at time t is found from (|29|) by 
taking squares, averaging over 6^0, Jq, dividing by 2|ep, and finally letting e — 0. With 
the help of (|E.71 ), this calculation yields: 

"^"^ sin^ (Jot/2) e 



{Et,o) = ^ / djQ . ; ; / = -t . (30) 
Sir sm (Jo/2) 4 
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The small contribution of the initial quasi-momentum in the atom's energy was 
neglected. Apart from that, (|30| ) is the same result as was found by the exact quantum 
mechanical calculation performed at e = and £ = 1 in section p.2.3| for the case of a 
uniform QM distribution, see (p3D. Thus the e-quasi-classical approximation reproduces 
the quantum behaviour at exact quantum resonance. 

The integral over Jq in (PD| ) collects contributions from all the invariant curves 
Jo =const. of the eSM at e = 0. Of these, the one at Jo = 27r leads to quadratic energy 
growth because it consists of (period 1) fixed points. This is called a classical resonance. 
It is responsible for the linear growth of energy (pO|), because the main contribution to 
the integral in (^) comes from a small interval ~ 27r/t of actions around Jo = 27r. Note 
that Jo = 2tt corresponds to /5o = 1/2, the "resonant" value of quasi-momentum. It is 
hence seen that the e— quasi-classical approximation explains the quantum resonances 
of the KR in terms of the classical resonances of the Standard Map. 

We shall now estimate (Et^^) for |e| > 0. The |e| > dynamics is maximally 
distorted with respect to the e = one for Jo in the vicinity of 2nn, that is, in the 
very region which is mostly responsible for the linear growth of energy at e = 0. Being 
formed of period-1 fixed points, the Jo = 2n7r, e = invariant curves break at |e| > as 
described by the Poincare-Birkhoff theorem |^ . The motion is then strongly distorted 
inside regions of size (in action) 6Jres astride J = 2mT. Such regions are termed the 
"main resonances" of the eSM, and a well known estimate has 5Jres ~ 4(fc|e|)-'^/^ ||41|| . 



Inside such regions, the approximation ( p9D fails quite quickly, so their contribution 
{Et)res in the mean energy has to be differently estimated. 

In the remaining part of the e— classical phase space the motion mostly follows 
KAM invariant curves, slightly deformed with respect to the e = ones, still with the 
same rotation angles. The contribution (-E'f)KAM of such invariant curves in the mean 
energy is therefore roughly similar to that considered in the integral (|30|) , provided Jo 
is therein meant as the rotation angleQ. On such grounds, in order to roughly estimate 
(Et^e) we remove from the integral (|30|) the contribution of the resonant action interval 
near Jo = 2tt, and we replace it by {Et)res'- 

{Et,,) ~ - m + {Et)res , (31) 



where 



SJres/^ Gin ^ 



8^ J~sjres/2 sm (J72) 
and J' is the deviation from the resonant value 2n. The contribution {Et)res may 
be estimated by means of the well-known pendulum approximation ETI]. Near the 



+ Higher resonances appear near all values of Jo commensurate to 27r. At small |e| such higher-order 
resonances affect regions of phase space, that are negligibly small with respect to the main resonance. 
They are altogether ignored in the present discussion. Also note that structures that are small compared 
to the "Planck constant" |e| are irrelevant for the purposes of the e— classical approximation. 



QR 



14 



eSM resonance, the motion is described (in continuous time) by the following pendulum 
Hamiltonian in the canonical coordinates J',i): 

Hres = ^iJ'y+MkcOs{l^). (33) 

The resonance width SJres is estimated by the separation (in action) between the 
separatrices of the pendulum motion. The period of the small pendulum oscillations 
is 2'Ktres where tres = so we use tres as a characteristic time scale for 

the elliptic motion in the resonant zone. One may altogether remove |e| from the 
Hamilton equations, by scaling momentum and time by factors (A;|e|)^^/^ = A/SJres, 
(A;|e|)^/^ = 1/tres respectively. Therefore, 

= ((j; - J'.Y) ~ k\e\G{t^\) , (34) 

for an ensemble of orbits started inside the resonant zone, where G{.) is a parameter- 
free function, whose exphcit expression involves elliptic integrals. This function results 
of averaging over nonlinear pendulum motions with a continuum of different periods, 
so it saturates to a constant value when the argument is larger than ~ 1. At small 
values (<^ 1) of the argument it behaves quadratically. The contribution to the total 
energy is then obtained on multiplying (^) by |e|~^5Jres/(47r), because only a fraction 
~ SJres/ (27r) of the initial ensemble is trapped in the resonant zone. As a result 

{E,)res ~ n-'\e\-'/'e/'G{t^\) . (35) 

When 6Jres is small, sin^(J'/2) may be replaced by J'^/4 in the integrand in (p2D, leading 

to 

$(t) ^ _t $o(tVW) 

with 

2 r sin^(s) 
<l>o(a;) = - / ds 



Replacing in (|3T|), we obtain: 

R{t, k, e) ^ ~ Hix) ^ 1 - $o(x) + —G{x) , 

X = t^/k\e\ = t/tres ■ (36) 

Hence R(t,k,e) depends on t,k,e only through the scaling variable x = t/tres- The 
width in e of the resonant peak therefore scales like [kt"^)^^. The scaling law (^) is 
demonstrated by numerical data shown in Fig. ^. The function H{x) was numerically 
computed: in particular, G{x) was computed by a standard Runge-Kutta integration 
of the pendulum dynamics (^). The scaling function H{x) decays proportional to x~^ 
at large x, because so do 1 — $o and 4:G{x)/{nx); the latter, due to the saturation 
of G. From Fig. |^ it is seen that $0 is quite slowly varying at a; > 4. The 
structures observed in that region are then due to G{x), which describes the resonant 
island. Numerical computation shows that G{x) saturates via a chain of oscillations 
of decreasing amplitude around the asymptotic value. These give rise to three local 
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maxima in the graph of x~^G{x), followed by a chain of gentle oscillations in the tail. 
The first and most pronounced maximum lies in the small- a; region, and is not resolved 
by the scahng function H{x), apparently because it is effaced by the rapid decay of 
1 — ^o{x). The subsequent maximum and its symmetric partner at e < are instead 
resolved and precisely correspond to the "horns" the right-hand one of which is marked 
by the arrow in Fig. ^. The oscillations in the tail of x~^G{x) are also well reproduced 
in the tail of H{x). 

The scaling law shows that at given k, t the energy curve (Et^e) vs. e decays 
proportional to |e|^^/^ past the "horns". As t increases the horns rise higher while 
moving closer and closer to r = 27r , because they are located at a constant value of 
X = t^yk\e\. Thus an overall |e|~^/^ dependence eventually develops (cf., Figs. ||). 
In the case when the smooth initial momentum distribution includes values uq ^ 
and/or is appreciably non-uniform in quasi-momentum, the statistical weights of the 
various phase-space regions are different. Scaling in the single variable t/tres still holds, 
but the scaling functions $o and G may be different. 

This analysis shows that the structure of the resonant peak is essentially determined 
by the main resonant island of the eSM. It neglects higher-harmonics resonances of the 
eSM, higher order islands, and especially the growth of the stochastic layer surrounding 
the main resonance. Such structures grow with |e| and are expected to introduce 
deviations from the scaling law (0). Hence this analysis is valid only if |e| is not 
too close to the threshold for global chaos \e\cr ~ l/fc. On trespassing this threshold, 
the critical regime of the eSM is entered. No isolating KAM curve survives, so the 
energy curve rises in time for |e| > \e\cr- Estimating the mean energy at relatively 
short times and for \e\k < 4.5 is difficult, because of the coexistence of unbounded, non- 
homogeneous diffusion and of elliptic motion inside residual stable islands. The increase 
of the curve with |e| at constant t is due to the decreasing size of the latter islands, and 
to the rapid increase of the diffusion coefficient (proportional to (|e| — |e|cr)'^, 7 — 3 at 
large enough t E^,E3i). 



4.1.2. Validity of the e— classical approximation near resonance. The e— quasi-classical 
approximation is trivially exact at all times for e = 0, as shown above. At nonzero e, 
it is valid for not too large times t, and it is in the long run spoiled by quantum, non 
e-classical effects. At |e| < \e\cr the e-classical motion is bounded by KAM curves, 
so the main quantum mechanism leading to non-e-classical behaviour is tunnelling. 
Estimating the related time scales is a non-trivial problem, because the 27r-periodicity 
in action of the e-classical phase space may enhance tunnelling, and even result in 
delocalization, depending on the degree of commensuration between 27r and the "Planck 
constant" |e|. For instance, if |e|/27r is rational, then the quantum motion will be ballistic 
asymptotically in time. This is just the ordinary quantum resonance of the quantum 
kicked rotor (section §). In order that one such resonance with |e| = 4np/q exists at |e| 
less than some |e|o, it is necessary that g > 47r/|eo|. It will show up after a time roughly 
estimated by |e| times the inverse bandwidth. The bandwidth is estimated to decrease 
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faster than exponentially at large q |]^,|T0|! so one may infer that the time of validity of 
the e-quasi-classical approximation is at least exponentially increasing with l/|e| as the 
exact resonance at e = is approached. 

At |e| > \e\cr the e— classical motion is unbounded, and the difference between e- 
classical and quantal energy curves vs. r is basically set by various quantum localization 
effects, including localization by Cantori close to the \e\cr 0,0- As a consequence, 
if t is large enough, then the e-classical curve lies much higher than the quantum one. 
Nevertheless the latter still rises with |e| at constant t, because of the growth of the 
localization length (that goes along with the growth of the e-classical diffusion |^). 

5. Decoherence induced by spontaneous emission. 

Loosely speaking, decoherence induced by external noise is expected to drive a quantum 
system towards classical behaviour. One in particular expects the quantum resonances 
for (5— kicked particles to be strongly impaired, whenever the noise spoils the conservation 
of quasi-momentum, because that purely quantal conservation law plays a key role 
in quantum resonances. This is exactly the case when decoherence is induced by 
Spontaneous Emission (SE) |21|,|3^, 0, |y, A beam of nearly resonant light of 



wave vector may induce an internal transition in the atom. The absorbed photon 
is thereafter spontaneously re-emitted in a random direction, and the whole process 
results in a random change of the momentum of the centre-of-mass atomic motion. The 
statistical average of this change equals hkj'. Although SEs are already produced at a 
quite small rate by the standing wave which produces the kicks, they have also been 
experimentally introduced in a controlled way by switching on near- resonant, low- 

intensity laser beams immediately after each kick. The additional laser beams induced 
an intensity-dependent mean number ris^ of SE events per atom and per kicking period 
between and 0.2. During each kicking period r, the SE-inducing laser was active over 
a (physical) time Tse ~ 0.067r. In our analysis we neglect the delay between absorptions 
and subsequent SEs. We hence assume that the atom undergoes random momentum 
changes due to SE at a discrete sequence of random times. Here we are only interested 
in the projection of such momentum changes onto the kicking direction^. Both the 
SE times and the corresponding momentum changes are modelled by classical random 
variables, independent of the centre-of-mass motion of the atom. This assumption is 
reasonable as long relatively small; otherwise, the atoms may, in the average, 

be slowed down or cooled - a velocity-dependent effect, which cannot be accounted 
for under this assumption. On the other hand, at the high field intensity of a quasi- 
resonant laser beam which is required for large ng-E stimulated emission would prevail. 
The statistical atom dynamics may be studied by investigating the stochastic Hilbert- 
space evolution of the atomic state vector |^,|2^|. More specifically, we assume an 

* The effect of randomly induced deviations from the straight hne is, however, an interesting question 
that should be studied separately. 

This approach was first implemented for the study of randomized (5-kicked rotor dynamics in ref . pq] . 
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initial incoherent mixture of plane waves with a distribution f{po). We further assume 
that all the variables describing SE events are given and fixed. Their specification defines 
a realization of the random SE events. Then we compute the deterministic evolution 
of the atomic state vectors up to time t. The quantum probability distribution of an 
observable in the final state depends on the chosen realization and on the initial state 
as well; averaging over both yields the final statistical distribution for that observable. 
The time evolution of distributions and related averages is the object of the following 
analysis. 

5.1. Kicked dynamics in the presence of SE. 

Let {ipit)) be the state vector of the atom immediately after the t—th kick. Let the 
integer z/^ denote the number of SEs during the subsequent SE-inducing window. Such 
events are assumed instantaneous. Denote the delay (in physical time) of the j— th 
event with respect to the (j — 1)— th one, and dj the momentum change of the atom 
induced by the j— th SE. For notational convenience a 0— th fictitious SE is assumed 
to occur immediately after the t—th kick, and a {ut + l)-th one immediately before the 
{t + l)-th kick, with do = d^^^i = 0. 

The state vector immediately after the {t + 1)— th kick is, apart from an inessential 
phase factor, 

\^{t + 1)) = 5(5t)e"''''°'^^^e-*"(^+^*)'/2|^(t)) , (37) 

where: 

j=l k=0 j=0 

Note that the 1st and the 2nd operator on the rhs of ( |57| ) commute. The conservation 
of quasi- momentum is broken by the operator S{6t); however, it may be restored by 
means of the substitution: 

\m) = S{5o + 5i + ... + 5t.,)\m) (39) 

which is a time-dependent, momentum shifting gauge transform; the resulting gauge 
will be termed the stochastic gauge in the following. Replacing ( ^91) in (|37|) , and using 

one easily finds: 

\^{t + 1)) = e-''='=°^(^)e-'"(^+^'+^o'''5»)V2|^(i)) . (40) 

This evolution does preserve quasi-momentum, so reduction to /3— rotor dynamics may 
be performed as described in section ^ leading to 

\Mt + 1)) = Mt)\Mt)) ; Mt) = e-'^-^We— , (41) 
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where 

t-i 

Vt = P + Xt + Y.^'^- (42) 

In this way the stochastic evolution ( pSj) has been separated in two parts. One of these 
is described in ( pOf ) by the operator S; the other is the evolution (^T]) in the stochastic 
gauge. The former is just a translation (in momentum) by the total momentum imparted 
by SE during the considered time; this part is totally classical, so it is the latter part 
that encodes the coherent stochastic evolution. 

5.2. A theoretical model for randomized kicked rotor dynamics. 

In experiments, SE occurs at random times within SE-inducing time windows, and there 
is one such window immediately after each kick. We shall say that a "S'i? evenf occurs 
at the integer time t, whenever at least one SE occurs in the window following the t— th 
kick. In this way, the number of SEs may be larger than the number of SE events. The 
probability of a SE event is denoted psE- We then assume: 

(51) SE events occurring at different times are statistically independent. Hence, the 
random variables 5t (t = 0, 1,2, ...) specifying the total (projected) momentum change 
produced by SE in the t— th kicking period (eq. (|38|) ) are independent, identically 
distributed random variables. 

The following two assumptions are at once the simplest and the strongest possible 
ones. They were chosen because they greatly simplify the otherwise still possible, yet 
cumbersome exact solution. Their validity will be discussed in section |5.2.5| , where it 
will be demonstrated that the main results derived below remain unchanged under less 
stringent and more realistic assumptions. 

(52) The finite duration of the SE-operating windows is negligibly small compared to 
the kicking period. Hence, Sj^^ ~ r while sj ~ for j < ut, so that xt ~ in (|38|) , 
and Vt - Po + T/o^t in (H. Different SEs may occur in the same kicking period, 
each separately contributing to the total momentum change St recorded in that period; 
nevertheless, their separation in time is neglected. 

(53) The occurrence of a SE event results in total randomization of the quasi-momentum. 
Given assumption S2, this is equivalent to assuming that the conditional distribution of 
the variable 6t mod(l), given that a SE event occurs at time t, is uniform in [0, 1] (in 
our units). We further assume a zero mean for the distribution of 6t. While no further 
specification is needed for the formal elaborations below, in numerical simulations we 
shall in fact use a uniform conditional distribution in [—1/2, 1/2] (in our units). 

(54) As in the SE-free case, we assume the initial statistical ensemble to be an incoherent 
mixture of plane waves. 
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5.2.1. Random Walk in Hilbert space, at exact resonance. Assumption (S4) will not 
be used in this subsection. The results of this subsection may therefore be used under 
different choices of the initial ensemble. Let SE events occur after to = at integer 
times ti = Ao,t2 = + ^i? ■ ■ ■ ,tj = tj_i + Aj_i, ... The variables A^, (j > 0) are 
integer, independent random variables. Under assumption (SI) they are distributed on 
the positive integers n with probabilities: 

pin) = Pse(1 - PsbT'' forn > , p(0) = , (43) 

where ps-E is the probability that at least one SE takes place in one kicking period. For 
all integers t > we define A^^^ =max{j : tj < t}, the number of SE events occurring 
not later than time t, and A^^o = 0. The integer random variables Nt, t > define a 
Bernoulli proces^ffj. 

After such preliminaries, we set out to study the evolution in the stochastic gauge, 
as defined by eqs. (0). The quasi- momentum /5 of a /5— rotor is constant in time, and 
for each rotor rjQ = f3 (SE events are allowed immediately after kicks, and no kick occurs 
at t = 0). In the stochastic gauge, the random propagator from time to time t for the 
/?— rotor is given by the ordered product 

t~i 

The subscript 5* on the Ihs refers to the stochastic gauge. The one-step propagators 
Ui3{t) are defined in eq. (|l|). Similar to what was done in section at exact resonance 
T = 2TTi one may write (cf. eq. (|^)): 

U^(t) = e-*'=™«(^~)e-^«*^ = e-^'^=°'(^")7^(^,) (45) 

where C,t = '^d-i'^flt =t 1). Although rjt is not restricted in the interval [0, 1] in eq. (^), 
the resonance condition allows for to be taken in [— vr, vr] in (^5)). Under assumption 
(S2), has a constant value in between t = tj and t = tj+i, such that 

j 

= Co + 2vr£ ^ Sm mod(27r) , -tt < < vr , (46) 

m=0 

where 6m is the total momentum imparted by the m— th SE event. Hence, = ^Nt 
in between the t—th kick and the (t + 1)— th one. A realization of the SE events is 
assigned by specifying the values of all the SE random variables just defined, which we 
collectively denote by the shorthand notations 6 for the random momentum shifts and A 
for the random times of SE events. Once the final observation time t and the realization 
are fixed, for notational convenience we re-define Ai^^_-^ = t — tN^_j^. Replacing (|46| ) in 
(1^), and then in (0), we get 

4-1 

ft This process should not be confused with the continuous Poisson process that may be used for the 



SEs occurring within one and the same kicking period; see subsection 5.2.5 
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By repeated use of 

eq. (PI) may be rewritten in the form: 

^sAt) = Mo + U + . . . + U_ Je-^'=^(^'^'^'*) , (47) 

where: 



t-i 



F{e, S, A, t) = ^ cos(^ + J2 U)- (48) 

s=0 r=0 

We next define •jj = J^m^Q A^fm- Replacing s in (|48|) by s = j + / with j = Ng, and 
summing over j, / separately, 

Nt-i Aj-1 

F(^, I, A, t) = J] ^ cos(^ + 7, + IQ . (49) 

i=0 1=0 

We further define: 

Zj = ^ exp(z7j + irQ ; Zm = ^Zj ; Wt = Zn,_, , (50) 

r=0 j=0 

SO finally 

F{e, 6, A,t) = \Wt\ cos{e + aigiWt)) . 

Note that Wt here differs from Wt in eq. ([T0|). Let the initial state of the atom be a 
plane wave of momentum Po = Uq + (3q. Given a realization (5, A), we operate on the 
corresponding rotor state (^ with the propagator (47). The (random) state of the rotor 



at time t is given, in the momentum representation, by: 

Jo 

cpt = (n - no)aYg(Wt) - n^Nit-i) ■ (51) 

The distribution of momenta n at time t is 

P(n,t|no,/3o,5, A,iV) = Jl_^S^\Wt\). (52) 

This is formally identical to (|T^, but now Wt depends on the initial quasi- momentum 
(3q and on the realization of the SE events as well. It is a stochastic process and the 
random state (|5l[) performs a random walk in the rotor's Hilbert space. 

5.2.2. Statistical Averages. Computation of statistical averages requires averaging over 
the SE random variables (5, A) and over the initial momentum Po = ""-o + /^o- Under 
our assumptions all such variables are classical random variables. Expectations (resp., 
conditional expectations) obtained by averaging over such classical variables will be 
denoted £{.} (resp., £^{.|.}). For instance, £^{.|po} stands for the average over the SE 
variables alone, given the value of po (or equivalently of no, /3o, ^o)- 

The large— t behaviour of the stochastic process Wt that drives the stochastic rotor 
evolution is ruled by the large m behaviour of the process Zm- The properties of the 
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latter process are completely determined by the assumptions (S2),(S3). Together with 
eq. P6| ) assumption (S3) entails that the are mutually independent random variables, 
uniformly distributed in (— 7r,7r) (with the possible exception of ^o, whose distribution 
is defined by the initial ensemble). This fact has the following consequences, that are 
derived in [Appendix lj| : the complex variables zj are pairwise uncorrelated whenever 
j + k > 1; moreover Zj, are independent whenever \j — k\ > 2. Thus the process is 
a random walk in the complex plane, and the distribution of approaches an isotropic 
Gaussian distribution in the complex plane as m — oo. On account of the properties of 
the Bernoulli process, the process Wt at large t has quite similar features. 
Moments of Wt may be explicitly computed at all t. For instance 

Nt-i 

j,k=0 

Nt-1 

= J2 + MNt-i) Re{S{z,z;\A}) , (53) 

j=0 

where ([B.3|) was used, and x(-) is the characteristic function of the strictly positive 



integers. The variables Aj were defined such that Ylflo^ ~ ^! hence, with the help 
of (i3) and (iJ) we find: 

8{\Wt\'' I A} = A^(Ao) - Ao + t + 2x{Nt-i) ReU{A^) . (54) 

Now Ao = ti if ti < t and Aq = t if ti > t, with ti the time of the 1st SE event. 
Furthermore, Nt_i = is equivalent to ti > t, and ti is distributed according to (|i3|). 
Therefore averaging over the random SE times A yields: 

= t Prob{ti < t} + M{t) Prob{ti >t} + C(t,psE) 

= t{l - + 4^Mit) + C(t,psE) (55) 

where Qse = 1 — Pse, and 

t 

C{t, Pse) = J2 P(^) [-^ in)-n + 2ReAr{n)] . (56) 

n=l 

With a smooth initial QM distribution (^) it is easily seen from the definitions 
and ( |B.4|) of M and Af that 



|C(t, Pse) I <3(5/o) 5^ np(n) , (57) 

n=l 

where (Sfo) is the maximum of |/o(/5) — 1| in [0, 1). In the case of a uniform initial QM 
distribution, Ai(t) = t and C{t,psE) = 0, so 

S{m'} = t, (58) 

like in the case without SE. With a non-uniform initial QM distribution, letting p^E — > 
at fixed t causes the 1st and the 3d term on the rhs of ( |5^ ) to vanish. If instead t cxo at 
fixed psE > 0, then the 1st and the 2nd term on the rhs approach t and zero respectively, 
exponentially fast; the 3d term remains bounded according to (|^). So the result which 
is obtained with a uniform QM distribution is asymptotically approached. 



QR 



22 



5.2.3. Growth of energy. Assuming that the initial state and the SE reahzation are 
given, and denoting 5t = J2o~^ "^s, the quantum expectation of the energy of the atom 
at time t may be written as: 

m = l I dp p'\{p\m)\' = 1 1 dpip+~5tr\{p\m)\' 

= 11 dpp'\{p\m)\' + ls^ + St I dpp\{p\m)\\ (59) 

where (^9]) was used. This expression has to be averaged over the initial statistical 
ensemble and over all SE realizations. Then 

S{~6n = D{t-l) , S{~6t} = 0, (60) 

where D = S{6^} is the mean square momentum change per period due to Spontaneous 
Emission. For an initial plane wave (^ of momentum Po = + /^o? with the help of 
(^) one finds 

dpp\{p\m)\'= E f dp {n + P)\{n\Mt))\ 



|2 



= J2in + Po)JLnoikm\)=no + /3o. (61) 

n 

where Jn(.) = (— )"J_„(.) and J2n'^ni-) = 1 were used. Similarly, 

J dpp'\{p\m)\' = E(^ + ^o)V^„„(A:|^^,|) , (62) 



whence, using (|K 



dpp'\{p\m)\' = ^k^mi' + {no + Pof. (63) 

Replacing (^) in (|5^) , the expectation of the last term in (p9|) vanishes due to (|60|) . 
The expectation of (|63D is found with the help of (|55D . Thus finally: 



£{E{t) - E{0)} = h'[t{l - gL) +M{t)ql, + C{t,p,^)] + \D{t - 1).(64) 



This result reduces to the SE-free one eq. (p3|) for q^^ = 1. The term on the rhs 
which includes fc^ as a factor is the mean energy in the stochastic gauge. With a 
uniform QM distribution, it reduces to (cf. (|58|)). so it does not contain any SE- 

related parameters, and is in fact identical to the result obtained in the SE-free case 
(for an initially uniform QM distribution). In experiments, D is typically small (see 



section |5. 2. 5| ), and the QM distribution is practically uniform; so the growth of the mean 



energy is but weakly affected by SE. However, a similar, albeit cumbersome computation 
of higher-order moments would reveal sharp differences, which reflect totally different 
ways of spreading of the momentum distribution in the two cases. At large time, this 
can be analyzed in detail, as shown in the next section. 

Under assumption (S3), the QM distribution of an atom is immediately turned 
uniform by the first SE event. The time scale for uniformization of the QM distribution 
is then tc = —1/ log(gsE)- Eq. (|64|) shows that for t ':$> tc the growth of energy is linear 
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with the coefficient + D/2, like in the case of a uniform QM distribution. On the 
other hand, since C(t,psE) is bounded in time, for (5/o)tc <^tc the growth of energy 
is dominated by the term M.{t), which is the same as in the SE-free, non- uniform case. 

5.2.4- Asymptotic momentum distribution. In this section we assume no = 0; averages 
over initial distributions allowing for n^ ^ may be easily implemented on the final 
results. We denote P{p,t) the momentum distribution at time t. We show that, as 
t — > oo, P{p,t) approaches a Gaussian distribution with mean value 0, in the sense 
that, for an arbitrary smooth function (f){p), 

lim(0)t= limVi / dp P{pVt,t)(f){p) = / dp (j){p)gD+ky2{p) (65) 

t— >oo +00 / / 

where Q^2 (p) denotes the normal distribution with zero mean and variance cr^ . To this 
end we compute: 



(0)i = ^|y dp\{p\m)\'<Pip/Vi) 
= sy dp\{pmw<i>i{p+~St)/Vi) 

= £: 4{k\WtM{n + /3o + 6t)/Vi) \ . (66) 



For t ^ 1 one may neglect corrections of order 1/vt in the argument of the smooth 
function 0, so 

{<p),^ £ i^-^jiikmrnn+M/vi) 

where [.] denotes the integer part. Asymptotically as t — oo, the statistics of Wt are 
determined by the fractional parts of sums of many 6s (cf. eqs. (^6D,(|50|)). Such sums 
of a large number of independent terms have a broad distribution, so their integer and 
fractional parts tend to be independent of each other as the number of terms in the 
sums diverges. The squared Bessel functions and the function in the last equation 
may then be separately averaged. Denoting 

Mp)^£mp+MVi)} (67) 

we may write 

^J2^{Jnik\Wt\)}Mn/Vi). (68) 

n 

As t — >■ oo, the distribution of Wt approaches an isotropic Gaussian distribution in 
the complex plane centred at 0. The variance is found from (|55D to be ~ t at 
large t. Hence the distribution of p = iVFfl is asymptotically at large t given by 
dFt{p) = 2t~^pdp exp(— p^/t). Consequently, 

POO 

(0),^ / dFtip)J2jlikp)Mn/Vi) 
Jo 
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= 2 1 dx xe-'''y2-JnikxVt)(f)tin/Vt). (69) 

The integral over p is a classical expectation, but the sum over ri is a quantum 
expectation instead. It may be written as: 

I,{kx) = Jl{kxVt)ct>t{n/Vt) = (O|00t(ri/2_^)/C|O), (70) 

n 

where 

If we regard as the Planck constant, then t ^ cxd is equivalent to a classical 

limit. In that limit t"^^^// corresponds to (angular) momentum p, and )C corresponds 
to p ^ p + kxsm{6). Therefore, the "classical" limit {t — > oo) for the momentum 
distribution in the state /C |0) is given by the distribution of kx sin(6'), with 6 uniformly 
distributed in [0,27r]. Replacing the quantum expectation (^) by the average over the 
related classical distribution yields 

\\mlt[kx)= / — 0oo(«;a;sm6') = / (72) 



J-TT 271" J^kx yA?2^^j2T7^ 

Substituting this in (|69D, and computing the elementary integral over x gives 



t— >oo 

On the other hand, by eq. 



oo -| 

lim (0)^ = / dx xe-"'joo(A;x) = I dp 0oo(p)e-^'/'' . (73) 



k^/7T 



Ap) = j dp' <p{p - p')Qd{p') 



where Qd is the limit {t oo) normal distribution of 5t/\/t. Recalling (|7^) and 
the definition of {(j))t given in ( |53D we immediately obtain the result claimed there. 
Hence, P(p, t) is asymptotically equivalent to a Gaussian with zero mean and variance 
k'^t/2 + Dt. Being just the leading term in the asymptotic approximation as t — oo, this 
misses those terms in the exact result ( |64|) which are bounded in time. The way (^) 



was derived from (|68| ) shows that decoherence turns the dynamics classical by causing 



the effective Planck's constant to decrease with time. An exact derivation of (l73D from 



is given in [Appendix C 



5.2.5. Discussion of the model, and numerical results. We shall now discuss 
assumptions (Sl-3); regarding (S4) see section |3l2 . 



(SI): The experimental time window for SE: Tse = 0.067r = 0.424, for r = 27r (or 
4.5/isec, in the Oxford experiments) is on the one hand very large compared to the time 
scales given by the Rabi frequency of the used atomic transition and by the inverse SE 
damping rate P^,|3^. On the other hand, Tse is very small compared to the kicking 



period. It is therefore reasonable to neglect memory effects inherent to SE processes [14 
so this assumption appears legitimate. 
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(52) : For the case when SE is induced by the kicking wave itself, this assumption 
remains vahd as long as the 5— kick approximation is valid. For the case of the Oxford 
experiments it may be to some extent supported by the smallness of Tse/t ~ 0.067. 
It should however be mentioned that the main results of our analysis: ( |6^ and 
the asymptotically Gaussian distribution still hold in the absence of (S2), though a 
considerably more involved analysis is required. Rather than delving into such analysis, 
we support this claim by numerical results to be presented below. 

(53) : Complete randomization of quasi- momentum after each SE- inducing cycle occurs 
when the mean number of SEs per period is large, ^ 1: then effective averaging leads 
to quasi- independent rjt in (^21) . If else rigE is small, then complete randomization requires 
that the conditional distribution of 6t mod(l) (given that at least one SE occurs in the 
t—th kicking period) be uniform in [0, 1]. On the other hand St is the sum of a random 
number of momentum changes due to single SEs. If these are assumed independent 
and identically distributed, then each of them has to be uniformly distributed in some 
interval of integer width. For a single transition in a 3-dimensional atom, the probability 
distribution of momentum shifts produced by SEs is not isotropic |^ . This in particular 
implies that the distribution of single SE, projected momentum shifts 6p is not uniform. 
In the case when the SE-inducing beam is orthogonal to the kicking direction, it has 
the parabolic form: 

n(s,)J<i-'^) ' l^fl^^-/^*^ (74) 

, otherwise, 

where C is a normalization constant, and -L are the (assumed to be orthogonal) 
wave vectors of the SE-inducing light and of the kicking light, respectively. This 
distribution is derived for a situation where SE from a Am = ±1 atomic transition 
is induced by circularly polarized light |^9|. The allowed change in momentum 6p 
is restricted within the interval [—kT/2kL,kT/'2kL], with kx/^kL — 1/2 (resulting in 
C ^ 1) in |3^,|3^. With a non- uniform distribution such as (|7^ some correlation is 



established between QMs in different kicking periods. The mean momentum change 
due to absorption followed by SE is hkj- for a single SE-inducing beam with wave 
vector kx- Our assumption of zero mean (along the kicking direction) is justified either 
when kr ki, or when the experimental arrangement uses two or more appropriately 
directed beams, whereby the atoms may be excited with equal probability. In such 
cases, the distribution of the projected 6p is more complicated (and closer to uniformity) 
than (|71D. Since the experiments use a large ensemble of Ce atoms, and SEs involve 
several hyperfine sublevels |3^,|2^, the assumption of a nearly uniform distribution of 



momentum changes seems, however, most appropriate. We shall nonetheless use (0) as 
a term of comparison in "type (II) simulations" (see below) in order to test the effects 
of deviations from uniformity. Such numerical data demonstrate that our assumption 
of a uniform distribution in an exactly integer interval of allowed momentum changes 
does not affect the results, for experimentally relevant times at least. 
For the uniform distribution of 6p in the interval [—1/2,1/2], = 1/12, so the 
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coefficient D in is D = n^^/l2. With the distribution (V) = 3/40, and 
D = ngE3/40. In the theoretical model based on assumption (S2) the distribution of the 
random times at which single SEs occur within one kicking period is totally irrelevant, 
so and Uqe enter as independent parameters. They have to be related to each other 
in order to make contact with experiments. A seemingly natural way assumes a Poisson 
distribution for the SEs occurring within one operating window, at least for not too 
large psE- In that case, psE = 1 — exp(— ?t,se)- 

We have performed numerical simulations of two types. Type (I) used all assumptions 
(Sl-4); type (II) had assumptions (S2),(S3) replaced by more realistic ones, allowing 
e.g. for free evolution in between successive SEs occurring in the same kicking period, 
and using various distributions of 6p, such as (|7^ , as discussed above. Type (I) 
simulations serve as a demonstration of the theoretical exact results, and much more 
as a term of comparison with type (II) simulations. The essential agreement between 
the two types demonstrates that our theoretical conclusions remain valid, under less 
stringent assumptions. Both types of numerical results were obtained by independently 
evolving rotors in a given Gaussian ensemble, and by incoherently averaging the final 
results. Random SE events were simulated as follows. After choosing values for Tse 
and psE = 1 — exp(— rigE); random SE times were generated in each kicking period 
from a Poisson distribution with the characteristic time Tse/'^se within the time window 
{tT,tT + Tse). To each random time a random momentum jump was associated, from 
the chosen distribution (uniform or parabolic). In type (I) simulations, such jumps 
were added to the quasi-momentum the rotor had at (integer) time t. The integer part 
of the result determined a corresponding shift in the computational basis of angular 
momentum eigenstates. The fractional part was used as quasi-momentum for a full 
one-period free rotor evolution. In type (II), free evolution was allowed in between 
subsequent SE times. In all cases the computational basis of momentum eigenstates was 
chosen as large as possible in order to model as faithfully as possible the ideal models 
analyzed in previous sections. 

Fig. ^ shows a long-time plot for different rates psE = 0.05 . . . 0.2, and for the two 
cases: type (I) with SEs happening immediately after the kicks (a), and type (II) with 
SEs within a finite time window (tse = 0.067r = 0.424 |3^) (b). For = 0.2 data 
is given in Fig. |^ (c) also for the parabolic distribution ([7^ . The energy growth is 
in all cases hnear with the predicted slope Ddec — + D/2, as discussed above 
(see eq. (0)). Fig. ^ presents the coarse-grained momentum distributions P„(t) (see 
section p.2.2| ) defined as the probability that the momentum p of an atom at time t lies 
in [n,n + 1) (in our units). They are computed for r = 27r and different SE rates, SE 
events immediately after the kicks in type (I), within the time window t^e = 0.424 with 
a uniform distribution, or SEs immediately after kicks with distribution (0) in type 
(II). With added decoherence, the distribution keeps spreading as a whole all the time, 
looking more and more Gaussian-like while it flattens out. Further remarks are given in 
section ^ Apart from statistically induced fluctuations in the wings of the distributions 
no significant difference between the different simulations (Figs. ^0) is detectable, and 
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our conclusions from the preceding subsections remain valid in all cases. We conclude 
that the results obtained in the work are not very sensitive to assumptions (Sl-4) that 
made the analytical treatment possible. 

5.3. e— quasi- classical approximation, in the presence of SE. 

The e— quasi-classical approximation introduced in section^ for the study of the coherent 
nearly resonant quantum motion is easily adapted to the model in the presence of 
SE, because the effects of SE were modelled by a totally classical noise. In the 
stochastic gauge, the e— classical approximation may be implemented in the /3— rotor 
propagators ( plD much in the same way as in subsection The resulting e— classical 
map corresponding to (pT]) is: 

It+i = It + ~k sin(^t+i) , Ot+i =et±It + ni + Tr]t . (75) 

We now exploit assumption (S2) and write rjt = (3 + 6t, where 6t = X]l=i is the total 
momentum imparted by SE up to time t. In order to turn off the stochastic gauge, we 
need to recover the accumulated SE momentum change, hence we change variables to 
Jj* = It + |e|5t. The momentum of the atom at time t is then |e|~"'^/i* + (3, where (3 is the 
initial quasi-momentum. Denoting rj^ = rjt + e/5/(27r£), a straightforward computation 
yields: 

= i: + \e\St+i + ~ksm{9t+i) , 
9t+i =9t± i; + 7r£ + 2niri; , 
Vt+i = Vt + ^t+i , 

r,l =r/?/(27r£). (76) 

The 5t are independent random variables, whose distribution is determined by the 
statistics of Spontaneous Emission. Numerical simulations of such noisy e— classical 
maps are shown in Fig. |], and very well match with the quantum computations 
at small |e|. The structure of the resonant peak in the presence of SE may be 
analytically analyzed, using ideas developed in section |. Under the substitution 
Jt = ±/* + 7t£ + 271 irj^ the map (|76| ) reduces to a noisy eSM , which differs from the 
eSM by a random shift T6t of the action J at each step. We assume an initially uniform 
QM distribution. At any SE time tj, the distribution of the ensemble in the phase 
space of the eSM is reshuffled by the random action change. Under the assumption 
of homogeneous distribution of single SEs in an interval of integer length, the resulting 
distribution of Jmod(27r) is approximately homogeneous over the unit cell of the eSM. 
Such randomization may be assumed to wash out correlations between the past and the 
subsequent random dynamics. Hence the scaling (|36[) may be used to write the energy 
at time t as 

{Et,,) ~ ^($^ A, H{A,/tres)) + -D'n^^t , 
i=o 

where (.) stands for average over all the Bernoulli realizations of the times of SE events, 
is the average number of Spontaneous Emissions per period, and D' = nsE~^{St) i^ 
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the mean square momentum imparted by a single SE. If is sufficiently large compared 
to 1, one may replace the Bernoulli process by the continuous time Poisson process with 
the characteristic time tc = — l/(log(l — Pse)) = I/'^-se- This process has the delays A 
distributed with the density exp(— A/tc)- Its statistics reduces to that of the unit 
Poisson process (tc = 1) by just rescaling all times by the factor 1/tc- This entails 

Nt-i 

{ J2 A, H{A,/tres)) ^ ^tcQit/t„ tJUes) , (77) 
j=0 



where 



1 " 

Q{u,v)^-{Y,A]H{A]v)). (78) 

j=0 

The superscript 1 specifies that the average is now over the realizations of the unit 
Poisson process: each realization has the continuous time interval [0, u] divided in 
subintervals Aj by a random number of Poisson events. We are hence led to the 
following scaling law: 

{E,,)r^D'± + k\Q{^,^) (79) 

or, equivalently, 

2{Et,,) - D't/t, 



2^2^ Qiu,v) , U = t/tc , V = tc/tres ■ (80) 

The scaling function Q{u,v) may be explicitely written in terms of the function H{x) 
by a routinely calculation reported in [Appendix L)| : 

4:Q{u,v) = uH{uv)e~''+ dx e-''xH{xv){2 + u - x) . (81) 

Jo 

Limiting behaviours of the scaling function Q{u,v) immediately follow from this 
equation, or from (^) itself. On one hand, for u = t/t^. ^ 1 the rhs in ( ffSD is a 
sum of a large number ~ t/tc of terms. In that limit, such terms are quite weakly 
correlated and may be independently averaged, leading to: 



oo 



M > 1 : Q{u,v) ~ / dx H{vx)xe ^ . (82) 

On the other hand, ioi t/tc -C 1, the sum reduces to the single term j = 0, with 
Aq = t/tc] hence 

M < 1 : Q{u,v) ~ ^uH{uv) . (83) 



In particular, ( p3D shows that ( [79|) coincides with (|36[) in the SE-free limit tc ^ oo. In 
the opposite limit, (^) shows that if k is fixed then the width in e of the resonant spike 
does not shrink any more with time when t ^ tc, and its width thereafter scales like 
{t1k)~^. The spike is therefore erased (that is, it is absorbed in the background) in the 
strong noise limit tc ~ 1. 

If fo{Po) is smooth but not uniform, the scaling in absence of SE of the form (|36D holds 
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but with a different scaling function H as explained in subsection [4.1.1| . Therefore the 
arguments of the present subsection leading to ( |79D should hold also in this case. 

Numerical simulations in Fig. ^ satisfactorily support this scaling law. Data are 
obtained in a similar manner as for the case without SE; however, one of the parameters 
u,v is varied, while keeping fixed either the other parameter or the ratio u/v. The 
theoretical scaling function Q{u, v) was numerically computed using in (^Tj) the function 
H{x) computed as described in section [4.1.1| . 



6. Reconciliation with experimental results. 

In the presence of SE induced decoherence, experimentally measured energies at fixed 
observation time tabs were found to exhibit resonant peaks near the resonant values 
T = 2ni (integer multiplies of the half- Talbot time |22|,^,|3^,^) that were higher than 
in the S13-free case. Numerical support for this observation was recently given in . 
Such observations may have been suggestive of an enhancement of quantum resonances 
due to decoherence: however, such a phenomenon has no match in the theory developed 
in the previous sections. This paradox will be resolved in this section. It will be shown 
that certain restrictions, that are unavoidably present in real experiments, depress the 
ideal resonant behaviour in a way, that is most severe in the absence of SE. So the 
explanation rather lies with the experimentally measured, SE-free peaks being lower 
with respect to the ideal case, than with the SE ones being higher. 

The most important experimental features not taken into account in the foregoing 
theoretical analysis are: 

(I) experimental kicks are not 5— like. The ideal model is then only valid as long as 
the distance travelled by the atoms over the finite duration tp of the kick is much 



smaller than the spatial period of the potential (tp = O.bfisec in \p% , |33|| ) . In our 
units Tp = tph{2kLy /M ~ 0.047. At large momenta this requirement is violated, 
and the atomic motion starts averaging over the potential. The small momentum 
regime is practically not affected by the replacement of the 5— function by a pulse 
of finite width. A proper theoretical description demands the 5— kicks in the atomic 



Hamiltonian to be replaced by appropriately shaped pulses [18|. For large momenta the 
pulses act adiabatically (if they are smooth), leading to trivial classical and quantum 
localization The classical phase space is then filled by KAM tori beyond 

some large momentum value Uref- For smooth pulses, the momentum Uref is inversely 
proportional to the duration of the pulse Tp, and the pre-factor of the proportionality 
depends on the shape of the pulse DTSf. The translational invariance in momentum 
required for quantum resonances is thus destroyed. The atom dynamics mimics the 
ideal lowest order resonances for a (possibly long) while |jl8|, but not the higher-order 
ones whose period (in momentum) is not very small compared to n^e/. This is an 
additional reason preventing experimental detection of high-order resonances, no matter 
how long the observation time. Fig. |l^ shows a simulation for an ensemble of rotors, 
with a rectangular pulse shape of width Tp. This does not include the smooth switching 
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on/off of pulses, as described e.g. in |2^; no substantial difference is however expected 
in the dynamics on relatively small time scales. In each kicking interval the rotors 
freely evolved over a physical time r — Tp. During the remaining time Tp they evolved 
according to the pendulum Hamiltonian (A/" + /?)^/2 + kcos{9). The latter evolution 
was computed by a Trotter-Kato discretization of the Floquet operator (equivalent to 
replacing the pulse by a thick sequence of 5— subkicks). 

(II) The experimental signal-to-noise ratio allows only a finite interval of momenta 
to be observed; in this border was ricut = 40 (data with counteracted gravity). 
Momenta with n > Ucut are not included in the experimental data of ||3^. Due to 
this fact, the theoretical momentum distributions have to be appropriately weighed 
prior to comparison with experimental ones. The crudest way is cutting the theoretical 
distributions beyond ricut and renormalizing the probability to 1. 

The effect of (I) and (II) on the ideal behaviour discussed in the previous section 
is easily understood in qualitative terms. We start with the SE-free case. The resonant 
growth of energy is stopped as soon as the ballistic peak in the tail approaches the closest 
of the two borders that are the effective cut-offs: (I) riref and (II) Ucut- If this happens 
earlier than the observation time, then the resonant peak is significantly depressed in 
comparison to the ideal case. We shall presently argue that such depression is mainly 
due to the cut-off (II) for the experimental cases of []33|,^. 

In the case of a rectangular pulse, n^e/ is not a precisely defined quantity, due to the 
slow decay of the Fourier harmonics of the pulse. It has to be meant in an effective sense. 
We hence resort to numerical simulations. In Fig. numerically computed momentum 
distributions are compared with those obtained in the ideal 5-kicked rotor case; according 
to such data, the effective riref should be located in the momentum range 70 — 120. 

The second cut-off Ucut is simulated by not counting momenta higher than ricut 
when calculating energies, momentum distributions etc. (the computational basis of 
momentum eigenstates is however much larger than Ucut)- Following experimental 
parameters we choose ricut = 40, the distributions are renormalized to unity 

after disregarding states with momenta larger /smaller than n = ±40. In Fig. |rT] the 
effect of this cut-off on the growth of the mean energy is shown. In the presence of Ucut 
the deviation from the ideal case appears somewhat earlier, as expected from ricut < nref] 
moreover, the deviation at t > 20 is strongly enhanced in the presence of ricut (and in 
the absence of SE). As shown in Fig. [l^ (a), the momentum distributions including 
both cut-offs (I) and (II) are stable in time, not moving at all in the centre around 
n = 0. The slight enhancement at \n\ ~ 15 — 40 as compared to the case without 
cut-offs (shown in Fig. |10D is only due to riref which to some extent acts like a reflecting 
boundary. The ballistic peak, however, which moves in momentum like n ~ 7rkt/2 
(see discussion after (|^)) is lost already after about t ~ 40/A; ~ 16 kicks, cf.. Fig. pi] . 
The peak is then beyond the cut-off (II). The estimated loss after about 16 kicks is in 
correspondence with the saturation of the mean energy vs. time at quantum resonance 
which has been observed in ||3^,|2^ for t > 15 in the experimental results, and in the 
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theoretical modelling as well (Fig. (4) in |p2|| ). 

The dependence of the mean energy on the kicking period r, which was shown 
in Fig. ^ for the ideal case of 5— kicks and no cut-offs, is strongly influenced by the 
cut-offs at exact resonance. This dependence in the absence of SE, with rectangular 
pulses and cutoff at ricut is shown in Fig. |13|(a). By comparing to Fig. |, we directly see 
that the only substantial difference is at resonant values r = 2tt,Att,6tt: cut-offs lead 
to lower resonant peaks. When the cut-off (II) ricut = 40 is applied in the ideal case 
of 5— kicks, no differences can be detected from the results plotted in Fig. |l^(a), again 
confirming that cut-off (II) is the crucial one. The resonance peaks are smaller, because 
the resonant growth of energy stops, as soon as the ballistically moving rotors hit the 
boundary (cf.. Fig. 0). Then the mean energy very quickly falls below its ideal value 



(after about 16 kicks in the plotted case), as can be seen in Fig. |TT 



Added SE totally changes this picture. The energy growth is now due to the overall 
broadening of the distribution, and not just to the ballistic peaks in the tail, as can be 
seen comparing the various parts of Fig. |^, where kicks are 5— like, and no cutoff is 
present. The distributions with weak SE are broader in the tails as compared to those 
with strong SE; the latter are however flatter in the centre, which is why they have 
roughly the same rms deviation. As already commented, in the S13-free case the quasi- 
momentum is constant in time, and atoms with quasi-momenta close to 1/2 travel faster, 
thus producing the long tails and the ballistic peaks at their edges. In the presence of SE, 
no atom may persist a long time in the fast-travelling quasi-momentum range, whence it 
is removed the sooner, the larger psE- Due to this reason, with SE the cut-offs are "felt" 
much later by the evolving distribution (Fig. [12], |^). Whereas the cutoff still prevent 
observation of the fastest atoms, they do not significantly affect the growth of energy 
until large times. Even then, the momentum distribution normalized within \n\ < ricut 
approaches the flat distribution in |?t,| < ricut (Fig. [T2|(b,c,d)), which has a limit value 
for the second moment significantly higher than the SE-free steady state distribution in 
the presence of the cutoff. 

In contrast to the SE-free case, the dependence of the mean energy on the kicking 
period r after 30 kicks is but slightly affected by the cutoffs when SE is present. This 
is shown in Fig. |l^ (b)(c), to be compared to Fig. ^b). In the experiments (Fig. (2) 



m 



32| , Fig. (6) in [Q), the peaks for all cases (a)-(c) are still smaller than in our Fig. 0, 
which can be explained by the extreme sensitivity of the resonance to 

all sort of perturbations besides those included in our present analysis, and also by 
difficulties in experimentally tuning to the exactly resonant values of r. Additional 
experimental restrictions, e.g., the experienced fluctuations of the potential depth and 
the resulting averaging over slightly different experimental realizations p3|,p7|,|i6[], may 
lead to a further reduction of the peak, especially in the case without decoherence which 
is most sensitive to any kind of disturbance. 
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Appendix A. Analysis of the Steady State Distribution (pO|). 

Appendix A.l. Proof of estimate 
Prom dH) it follows that 

1 r'^l'^ 

M*(t) = — ^ dx da Jni^) ) -2 := A; sin(a;) csc(a).(A.l) 

|n|>7V |n|>iV 

for any positive integer A^. In [Appendix A.2| we show that: 

|n|>JV 

We choose < e < tt and we use ( [A. 21 ) to bound the sum in the inner integral in 
( |A.1| ) when |a| > e/2; otherwise we use the upper bound 1. Noting that \z\ < 7rA;/(2e) 
whenever tt/2 > \a\ > e, 

where k' := /c7r/2. We now minimize the rhs by choosing 
ek' /87rA^2\^ 



2N V et ' (*-^' 
(which is indeed not larger than vr whenever N > kx 1.03 . . .). Replacing ( [A. 41 ) in ( |A.3| ) 
yields the estimate (PI|). 



Appendix A.2. Proof of ineq. I\A.2{) 



Using the bound ( [E.6|) and the power series expansion of Bessel functions |^ 

2n 



(JXJ (JXJ X 

n— n Ti— n \ 



z\e 



r/2 



n=0 n=0 

for any real r. It follows that: 



(A.5) 



n=N 



Ineq. ( [A.2|) follows upon optimizing with respect to r. 
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Appendix A. 3. Proof of the asymptotic formula (p^j. 

For z G [—1, 1] we define: 

f{z) ■= ^ r'\aJl{kzcsc{a)) ioi z ^ Q ; /(O) = i (A.6) 

The integrand in ( [A.6|) is meant = for a = 0, z 7^ 0. Using tfie integral identity (Pi).3| ) 
for Bessel functions, (pOD may be rewritten as: 
1 /•''^ 

M* = - dx cos(2na;)/(sin(x)) , (A.7) 

so, for \n\ > 0, M* is the 2n— th coefficient in the cosine expansion of /(sin(a;)). The 
function f{z) is different iable in [—1, 1] \ {0}. It will be presently shown that 



4k 

f{0+)= hm /(z) = ^ . 

z-*0+ TT 



Since f{z) is an even function, it will follow that its first derivative is discontinuous at 
z = 0. We choose e > and write 

f{z) = f,{z)+g,{z) , f,{z) := ^ [ da J',{kz csc{a)) . (A.8) 

Then is differentiable around 0, with g'^{0) = 0. Hence, /'(0+) = /^'(0+). Next we 
note that if z > 0, 

f^{z) = n'^z'^ J da F{kz csc{a)) , 

where F{x) := xJo{x)Jq{x). Noting that 

I csc(a) — a^'^l < c\a 
for < I a I < 7r/2 and some numerical constant Ci, one easily finds 

fXz) = -R-^z-^ j da F{kza-^) + 0(e) 

poo 

= -2k-n'^ j duu-^Jo{u)Ji{u) + 0(e), (A.9) 

where Jq{z) = —Ji{z) was used. Letting z — > 0+ and thereafter e 0+ we obtain: 

4k 

/'(0+) = lim /;(0+) = -2k7i-^ / du m-Vo(m) Ji(m) = . 

The integral was computed by using ( |li).5| ) and then formula 11.4.36 in |^ . 
Next we recall from ( [A.?] ) 

^ 00 
f{smx) = -M* + ^M*cos{2nx) . 

^ n=l 

According to the above analysis, the derivative of this function jumps by — Sfc/vr^ at x = 
jvr {j any integer). Hence the 2nd derivative has the singular part —Sk-K^"^ 5(x — jvr), 
leading to the asymptotic value — IGfcvr"^ for the coefficients in its cosine expansion. This 
yields 

^ 4k 

M„ ~ „ „ as n ^ 00 . (A. 10) 
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Appendix B. Statistics of the process Z^- 

Appendix B.l. Independence of the variables Zj. 

We show that, for any integers n, m, {m > n), the variables {zn, Zm) are independent 
of the variables (2:0, z^) whenever k < n — 2. . It suffices to show that {zn, Zm) are 
independent of (^o; •••;^fc,^0; ■■■^k)- To see this, let / be an arbitrary (Borel) function 
of m — n + 1 complex variables, and consider: 

Mk =£{f{Zn,...,Zm)\io,--,ik,Ao,-;Ak}- (B.l) 

Looking at (|50|), and recalling that the C,j are mutually independent, one notes that 
(zn, Zm) depend on Aj (0 < j < k) through the factor exp(i Ylo^j^j)^ hence only 
through Ylo^j'^j niod(27r). Therefore, ( |B.1|) is a function of the variable /i^ = |j Aj 
mod(27r) alone: Aik = A^fc(/ifc). Furthermore, since k + 1 < n, 



dP{^k+i, Afc+i) Mk+iifik + ffc+iAfc+i) (B.2) 



because C,k+i,Ak+i are independent of past variables; here dP{.,.) is their joint 
distribution. Now ^^+1 is independent of the integer A^+i, and it is uniformly distributed 
in (— 7r,7r). Then the integral does not depend on fik, so 

Zn, Zm)\^0, ^fc, Ao, Afc} — S{f{ Zyi , . . . , Zm ) } ) 

which proves the announced independence property. 

Appendix B.2. Correlations of variables Zj. 

We prove that, for any j,k > such that j + k > 1 

S{zjzl\A} = 6jkAj , S{z,Zk\A} = 0. (B.3) 



whereas 



S{z,z*\A}=^iAo) = £ / dPi^o)e''^' , 

j = l ^ 
J _7i- 

£{\zom = MiAo) ^ /ciP(eo) ""'?f/°9f (B.4) 

J sm (^0/2) 



where 



i-i . . 

rfP(^o) = (27r£)-Meo 5^ /o(/5,) , = 1^ + 2 + 7 • (^-^^ 

is the distribution of ^0 and /o is the probability density of the initial quasi-momentum 
(II)- 
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In order to show ( [B.3D , we denote: 



r=0 



SO that: 

j-l k-l 



1=0 m=0 

Let j > k, j + k > 1, Then j > 1, and the 1st product has the factor Hence ( p.6|) 

depends on via this factor alone (if j' 7^ + 1) or via this factor multiphed by <^*j-i 
(if J = k + 1), leading to a factor e*^^-i(^J-i~') with I < Aj_i — 1. In both cases averaging 
over yields zero because ^j-i is uniformly distributed whenever j > 1. The case 
j < k is recovered by complex conjugation. The 2nd equality in ( pj.3D is straightforward. 
If j = k, then from (|50|) it follows that 



— 7r 

where dP{^j) is the distribution of ^j. For j > 0, dP{^j) = d$,j/{2n), and the integral 



is computed according to (|E.71 ). The 1st equation in (|B.4D results of a straightforward 



calculation using the definitions (|50D of 2:0 and zi. 



Appendix B. 3. Central Limit property. 

The properties of the process allow to conclude that its distribution is asymptotically 
Gaussian, thanks to known results about the Central Limit Theorem for weakly 
dependent sequences Isotropy of the hmit Gaussian distribution easily follows 

from computing i^{R-e^(^^ e~*^2;j)} (the mean square displacement along the direction 
6* in iV steps). Using ( [B.3| ), the result is independent of 9. 

Appendix C. Gaussian asymptotics for momentum distributions. 

Throughout this appendix, the Fourier transform of a function / is denoted by /. 



Notations are otherwise identical to those in section |5.2.4| . In particular. 



Replacing this in (|5PD, and using the Bessel function identity ( pi).4| ), 
(0)t = -^ / dFtip) [ du(i)tiu)Joi2kpsm{u/2Vi)) 



'2ti 
2 

/2^ 



/•OO /» 

y dzxe-^'J du ^tiu)Jo{2kxVism{u/2Vi)) {C.l] 
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which in the hmit t ^ oo yields 

2 roo 



'271 







dx xe / du 4>ooiu) Jo{kxu) 



Substitution of 



{u) = -= I dp <p^{p)e-^P- (C.2) 
V Ztt 



and of ( [E.5|) yields the t ^ oo limit: 



POO PKX 

2tt~^ / dx xe~^ / (ip 

Jo J-kx 



Interchanging integrals, we obtain the same result found in the text (equation (^)) by 
different means. 

Appendix D. Derivation of eq. ( ^) . 

Let Ap {j = 0, 1, ...) be real nonnegative, independent random variables exponentially 
distributed with ^^{Aj} = 1. For j a nonnegative integer denote Sj = J2k=o^k, and 
s_i = 0. For given m > let A^^ =max{j : sj < u}. We shall compute the expectation 
of the random variable 

ATI 

/„^5^/(Aj) + /(n-.^.), (D.l) 

3=0 

where f{x) is a given nonrandom function; the sum in eq. ([78|) is of this form, with 
f{x) = xH{xv). We write = YlT=o f^,v where 

fu, ^ Xiu - s,^^)[xiu - - A])/(Aj) + 

+ x(A] + s,„i - u)f{u - , (D.2) 

and x(-) is the unit step function. Then, denoting G{x) = dsf{s)e~'^, 

£{fujs,^i} = x(r)[G(r) + fir)e-^] , r = u- . (D.3) 

Therefore, 

oo 

£{fu} = Y.^{£{fuAsi-i}} = G{u) + /(n)e-" + 

+ J2 dPj{x)[G{u-x) + f{u-x)e''-''] , (D.4) 

where dPj{x) = dx e~^x^~^/{j — 1)! is the distribution of Sj_i for j > 0. Summing over 
j and replacing the definition of G{x) we finally obtain 

S{f^] = 2j rfxe-^7(x) + /(M)e-"+ / dx e"" f {x){u - x) . (D.5) 
Jo Jo 

Eq. ( ^11) in the text is obtained on substituting f{x) = xH{xv). 
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Appendix E. Some formulae used in the text. 

The following formulas involving Bessel functions were used in the text and previous 
appendices, these are taken from |^] or derived from formulas there (we give the 
corresponding numbers in []): 



1 



a+2TT 



— I ci^e*^^°^V"' = I'^Jniz) [9.1.21] (E.l) 
27r ./„ 



E Jnix) = [9.1.76] (E.2) 

n=—oo 

/ dx Jl{bsm{x)) = / c/x cos(2nx)Jo^(6sin(x)) [11.4.7/8] (E.3) 
Jo Jo 

J2j^{z)e'^' = Jo(2^sin(t/2)) [11.4.8] (E.4) 

n 

T / N 1 /""^ , cos(2;x) , , s 

Jo{z) = - / dx ^ ^ ' [9.1.18] (E.5) 



\Jniz)\ < 1 



Z 



\n\ 



The following integral was used in calculating moments of energy: 

"^^ , sin^rtx) ^ 

dx ^ ; = 2TTt . (E.7) 



sm {X} 
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Figure 1: Time averaged stationary momentum distribution P*{p) for a uniform 
distribution of initial momentum po = (3q in [0,1), obtained by plotting (|T^ vs. 
p = n + f3o- The solid line drawn through the peak tops corresponds to a decay oc 1/n. 
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Figure 2: Evolution of coarse-grained momentum distributions (|T8D at resonance r = 27r, 
k = O.Svr for an ensemble of 10^ atoms, without decoherence. The initial momentum 
distribution is centred Gaussian, with rms deviation a ~ 2.7. (a): distributions for 
t = 30 (solid lines), t = 50 (dotted), t = 100 (dashed), t = 200 (dashed-dotted). (b): 
doubly logarithmic plot of the distribution at t = 1000 (solid line), compared to the 
asymptotic formula 4/c/(7r^n^) (dashed-dotted line) (^). No cut-offs are used, so the 
distributions characterize the ideal behaviour of an ensemble of (5-kicked particles. 
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Figure 3: Mean energy after t = 30 kicks, vs. the kicking period r, for an ensemble of 
10^ 5— kicked atoms, without momentum cut-offs, with the same initial distribution as 
in Fig. ^ and k = O.Stt. (a) no decoherence, (b) added spontaneous emission with rate 
PsE = 0.2. Step-size in r: 6t ~ 0.03, avoiding simple rational numbers; same results 
were obtained for a high-resolution random grid in the kicking period r. 
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Figure 4: Magnification of Fig. ^(a) near tlie resonance r = 27r. Quantum data taken 
from Fig. |^(a) (circles) are compared witli the mean energies of an ensemble of 10^ 
e— classical atoms (solid line) with the same initial momentum distribution, evolving 
under the e— classical dynamics (pT]). The value of r corresponding to the small peak 
on the right of the resonant spike is marked by an arrow for reference to Fig. ^ 
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Figure 5: Demonstrating the scaling ( pq ) of the resonant peak, in a right neighbourhood 
of r = 2tt. Open circles correspond to different values of the parameters e, k, t, randomly 
generated in the ranges 1 < t < 200, 0.001 < e < 0.1, 0.1 < A; < 50 with the constraint 
0.01 < ke < 0.2. In each case an ensemble of 2 x 10^ e— classical rotors was used to 
numerically compute the scaled energy R{t,k,e) (^61), with a uniform distribution of 
initial momenta in [0, 1] and a uniform distribution of initial 6 in [0,27r). Full squares 
present quantum data for k = O.Svr, t = 50 and t = 200. The solid line through the data 
is the scaling function H{x) of (^) obtained by direct numerical computation of the 
functions ^o{x) and G{x). The dashed line represents the function $o(a^); the dashed- 
dotted line has slope —1 and emphasizes the decay described in the text. The arrow 
marks the value of the scaled detuning x which corresponds to the arrow in Fig. H. 
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Figure 6: Average energy vs. time t at exact resonance r = 27r for the same ensemble 
of atoms as in Fig. for k = O.Svr, in the presence of SE events simulated in different 
ways as described in the text: (a) a random number of SEs occur immediately after 
kicks, each causing a momentum change Sp uniformly distributed in [—1/2, 1/2] (type 
(1) simulation), (b) SE times are Poisson-distributed in a window = 0.067r, with 
free evolution in-between them; 6p is distributed as in (a), (c) SE times as in (a), but 6p 
has the parabolic distribution (|71D with ki/^kx — 0.476. Rates of spontaneous emission 
PsE = 0.05 (solid), j9sE = 0.1 (circles), psE = 0.2 (plusses). The theoretical prediction 
(see text) for the coefficient of linear growth Ddec is approximately 1.59, whereas the 
data lead to D^ec ~ 1-58 — 1.60 except in (b) for p^E = 0.05 where it takes the value 1.55 
(strong fluctuations). The insets zoom into the region close to t = 200. No momentum 
cutoffs are used. 
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Figure 7: Evolution in time of coarse-grained momentum distributions for the same 
initial ensemble as in Fig. and for k = 0.8, in the presence of SE. (a) = 0.1, (b) 
PsE = 0.2, (c) psE = 0.8, for t = 50 and t = 200. The SE events are simulated in different 
ways. Solid hues were computed like in (a) in the previous figure; circles and squares, 
like in (b); plusses and crosses, like in (c) there. 
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Figure 8: Analogue of Fig. with the same initial ensemble, for k = O.Svr, in the presence 
of SE. Results of full quantum calculations (circles) and of e— classical ones (solid lines) 
in the presence of SE are compared near the resonance r = 27r, for different times and 
different rates of SE: (a) psE = 0.1, t = 30, (b) psE = 0.1, t = 50, (c) psE = 0.2, t = 30, 
and (d) psE = 0.2, t = 50. The quantum simulation was type (I), while the e— classical 
simulations used the map (|76D . 
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Figure 9: Demonstrating the scaling law ([7D| ) in a right neighbourhood of r = 27i. 
In (a), (b) the quantity on the Ihs of eq. ( pOD is plotted vs. one of the parameters 
u = t/tc or V = tc/tres while keeping the other fixed: (a) v = 2, (h) u = 4. In (c) the 
ratio u/v = 4 is fixed. Open symbols correspond to different values of the parameters 
t,tc,k,e, randomly generated in the ranges 1 < t < 200, 5 < tc < 60, 0.001 < e < 0.1, 
0.1 < A; < 20, with the constraints 0.001 < ke < 0.2 and t^Vke = 2 in (a), t/tc = 4 
in (b), t = Atl\/ke in (c). In each case an ensemble of 2 x 10^ e— classical rotors was 
used, with a uniform distribution of initial momenta in [0, 1] and a uniform distribution 
of initial 6 in [0,27r)). The random momentum shifts at each step of the e— classical 
evolution (^) were generated from the uniform distribution in [—1/2, 1/2]. Full squares 
represent quantum data for k = O.Svr, and e = 0.01 in (a), e = 0.05 in (c), and t = 50 
and t = 100 in (b). The solid lines correspond to the theoretical formula ([81|). 
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Figure 10: Coarse-grained momentum distributions for the same initial ensemble as in 
Fig. ^ and without SE, for k = O.Svr and r = 27r. The ideal case of 5— kicks (solid line) 
is compared to the case of rectangular pulses (open circles). Times are t = 10 (a), 30 
(b), 50 (c), 200 (d). 
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Figure 11: Effect of finite pulse width and of momentum cutoff on the growth of the 
mean energy for the same initial ensemble as in Figs. ^ and |10|, for k = O.Stt and r = 2tt, 
without SE (a) and with = 0.2 (b). Solid lines are for the ideal 5— kicks and no 
momentum cutoff; dashed lines for rectangular pulses, no cutoff; dashed-dotted lines for 
rectangular pulses and momentum cutoff at ricut = 40. In (a) the energy is significantly 
depressed by the cutoff after t > 20, not so in (b). 
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Figure 12: Evolution of coarse-grained momentum distribution for the same ensemble 
as in Figs. ^ and |1^, for r = 27r and k = O.Stt, with rectangular pulses and momentum 
cutoff at Ucut = 40. (a) = 0, (b) psj, = 0.1, (c) = 0.2, (d) = 0.8, after t = 10 
(solid lines), t = 50 (dashed), t = 200 (dashed-dotted) 
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Figure 13: Mean energy as a function of the kicking period r for the same ensemble of 
SE and the same initial distribution of atoms as in Fig. after 30 kicks for k = O.Stt. 
(a) no decoherence = 0, (b) = 0.1, (c) PsE = 0.2. The width of the rectangular 
pulse is Tp = 0.047 and ricut = 40, as in the experiments [32,3^. The shown range of r 
corresponds in laboratory units to 21.2 . . . 254.7 fisec. 



